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ABSTRACT 

The Spitzer /Infrared Spectrograph spectra of three spectroscopically anoma¬ 
lous galaxies (IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948) 
are modeled in terms of a mixture of warm and cold silicate dust, and warm and 
cold carbon dust. Their unique infrared (IR) emission spectra are characterized 
by a steep ~ 5-8 /irn emission continuum, strong emission bands from polycyclic 
aromatic hydrocarbon (PAH) molecules, and prominent silicate emission. The 
steep ~ 5-8 /irn emission continuum and strong PAH emission features suggest 
the dominance of starbursts, while the silicate emission is indicative of significant 
heating from active galactic nuclei (AGNs). With warm and cold silicate dust of 
various compositions (“astronomical silicate,” amorphous olivine, or amorphous 
pyroxene) combined with warm and cold carbon dust (amorphous carbon, or 
graphite), we are able to closely reproduce the observed IR emission of these 
galaxies. We find that the dust temperature is the primary cause in regulating 
the steep ~5-8 /irn continuum and silicate emission, insensitive to the exact sil¬ 
icate or carbon dust mineralogy and grain size a as long as a < 1 /irn. More 
specifically, the temperature of the ~ 5-8 /irn continuum emitter (which is es¬ 
sentially carbon dust) of these galaxies is ~250-400 K, much lower than that of 
typical quasars which is ~640 K. Moreover, it appears that larger dust grains are 
preferred in quasars. The lower dust temperature and smaller grain sizes inferred 
for these three galaxies compared with that of quasars could be due to the fact 
that they may harbor a young/weak AGN which is not maturely developed yet. 
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Subject headings: dust, extinction — galaxies: active — galaxies: individual 
(IRAS F10398+1455, IRAS F21013-0739, SDSS J0808+3948) — infrared: galax¬ 
ies 


1. Introduction 

Dust is a ubiquitous feature of the cosmos. It is present in a wide variety of astrophys- 
ical systems, including active galactic nuclei (AGNs) and starburst galaxies. Amorphous 
silicates and some form of carbonaceous materials (e.g., graphite, amorphous carbon, and 
hydrogenated amorphous carbon) are thought to be the two dominant cosmic dust species. 
Silicate dust reveals its presence in AGN and starbursts through the 9.7 and 18 /mi spectral 
features which respectively arise from the Si-0 stretching and O-Si-O bending vibrational 
modes (Henning 2010). The detection of the 3.4 /im C-H stretching absorption feature in 
AGN indicates the presence of aliphatic, chain-like hydrocarbon dust in AGN (e.g., see 
Imanishi et al. 1997, Mason et al. 2004). Starbursts lack the 3.4 //m absorption feature. This 
may be related to the fact that starbursts, particularly in their nuclear regions, are often 
heavily obscured by dust and ice. In the Milky Way, the 3.4 /im absorption feature is seen 
in the diffuse interstellar medium (ISM), but not in dense clouds where the 3.1 /im H 2 0 ice 
absorption feature is strong (see Pendleton & Allamandola 2002). Starbursts often exhibit 
strong H 2 0 ice absorption at 3.1 and 6.0 /im (e.g., see Spoon et al. 2004), while these ice 
features are not seen in AGN where it is too hot for ice to survive against sublimation (see 
Li 2007). Furthermore, starbursts also emit prominently at 3.3, 6.2, 7.7, 8.6 and 11.3/mr 
which are often collectively ascribed to polycyclic aromatic hydrocarbon (PAH) molecules 
(Leger & Puget 1984, Allamandola et al. 1985), however, these PAH emission bands are not 
detected in AGN (see Li 2007). The absence of PAHs in AGN is commonly attributed to 
the destruction of PAHs by extreme UV and soft X-ray photons in AGN (Roche et al. 1991; 
Voit 1991, 1992; Siebenmorgen et al. 2004). 

In starbursts, the 9.7 and 18 /im silicate features are always seen in absorption. In AGN, 
in contrast, the silicate features can be seen either in absorption or in emission , or can not be 
made out at all in the continuum. AGN are thought to be surrounded by an optically thick 
dust torus which would block the photons from the broad line region and accretion disk if 
they are viewed edge-on, through the torus (“type 2”), or would allow the detection of broad 
emission lines if they are viewed face-on (“type 1”). The anisotropic torus is invoked by the 
unification theory of AGN to account for the observational dichotomy (Antonucci 1993; Urry 
& Padovani 1995). In type 1 AGN, silicates are expected to emit at 9.7 and 18/mr where 
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hot dust can be directly detected, while in type 2 AGN, silicates are expected to be seen in 
absorption because of torus obscuration (see Li 2007, Mason 2015) @ Moreover, starbursts 
show a steeply rising emission continuum at ~ 5-8 /an before the onset of the 9.7 pa n silicate 
feature, while the ~ 5-8 pan. emission continuum is flat in AGN. 

Very recently, Xie et al. (2014) studied the ~ 5-40 pan infrared (IR) spectra of IRAS F10398+1455, 
IRAS F21013-0739, and SDSS J0808+3948 obtained with the Infrared Spectrograph (IRS) 
on board the Spitzer Space Telescope (Houck et al. 2004). They found that the IR spectra of 
these galaxies are anomalous. Spectroscopically, on one hand they resemble that of AGN in 
the sense that the silicate features in these galaxies are seen in emission ; however, they also 
exhibit strong PAH emission features which are absent in AGN. On the other hand, they are 
like starbursts in the sense that the ~ 5-8 pan emission continua of both these galaxies and 
starbursts all steeply rise with wavelength A and they all show strong PAH emission features. 

However, the silicate features seen in emission in these galaxies are often seen in absorption 
in starbursts. Furthermore, the steep ~ 5-8 pan emission continuum seen in these galaxies 
is much flatter in AGN (see Figure 1 of Xie et al. 2014). Let d In F u /d In A be the slope of 
the ~ 5-8 /ini emission continuum, where F v is the observed flux at frequency u, X — c/v 
is wavelength, and c is the speed of light. On average, d\nF v /d\n A ~ 0.8 for AGN, while 
for IRAS F10398+1455, IRAS F21013-0739, and SDSS J0808+3948, dlnF u /d\n A«4.1, 4.2, 
and 4.6, respectively. A detailed description of the spectral properties and a comparison of 
the Spitzer/ IRS spectra of these three galaxies with that of typical starubursts and quasars 
can be found in Xie et al. (2014). 

The focus of this paper is to explore the unique properties of the dust in terms of 
composition, size and temperature in these three spectroscopically anomalous galaxies. This 
paper is organized as follows. We describe the data in §[2j In j|3] we describe the dust 
model and calculate the mass absorption coefficients of various silicate and carbon dust 
species in the IR. The model-fitting to the observed IR emission is presented in §HJ In §[5]we 
compare the PAHFIT approach with the spline approach, and also compare the IR emission 
of our three galaxies with that of quasars and IRAS FSC10214+4724, an ultraluminous IR 
galaxy (ULIRG). The major conclusions are summarized in jj6j Throughout the paper, we 
assume a cosmological model with H 0 = 70 kms -1 Mpc -1 , Q m = 0.3 and fl A = 0.7. 

The L 0 represents solar luminosity of 3.826 x 10 33 ergs s' 1 and M 0 represents solar mass of 
1.989 x 10 33 g. 


Where are some exceptions: some type 2 AGN exhibit silicate features in emission instead of absorption 
(Sturm et al. 2006; Mason et al. 2009; Nikutta et al. 2009). 
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2. Observations and Data 

These three galaxies were found when cross-matching the Sloan Digital Sky Survey 
(SDSS) and Spitzer/ IRS low resolution spectra (L. Hao et al. 2015, in preparation). The 
Spitzer /IRS spectra of these galaxies are taken from the Cornell AtlaS of Spitzer/IRS Sources 
(CASSIS) which includes ~ 13,000 low resolution spectra of >11,000 distinct sources observed 
in the standard staring mode and provides publishable quality spectra (Leboutciller et al. 
2011; L. Hao et al. 2015, in preparation). We tabulate their basic properties in Table [0 

While the ~5-8/im continuum and the 9.7 and 18 /im silicate emission features result 
from dust grains which, heated by the AGN central engine, attain an equilibrium temperature 
(see Li 2007), the PAH emission features arise from PAH molecules which are stochastically 
heated by single stellar photons from their host galaxies (see Draine & Li 2001). Due to 
the differences in their emission mechanisms, carriers, and emitting regions, to facilitate a 
detailed examination of the dust emission, we first remove the PAH emission features and 
the ionic emission lines from the Spitzer/ IRS spectra of these galaxies. This requires an 
estimation of the dust emission continuum underneath the PAH and ionic emission lines. 
We take two approaches in defining the continuum. First, we approximate the continuum as 
a sum of starlight and blackbodies of different temperatures using a modified version of the 
PAHFIT software (Smith et al. 2007). The continuum determined in this way is hereafter 
referred to as “the PAHFIT continuum.” We also define a “spline continuum” which is 
obtained by selecting anchor points at 5-7, 14.5-15.0 and 29-30 /mi of each Spitzer/IRS 
spectrum to form a underlying continuum and then fitting these points with a spline function 
(L. Hao et al. 2015 in preparation). In both approaches, the PAH emission features are fitted 
with a number of Drude profiles which are expected for classical damped harmonic oscillators 
(see Li 2009), while the ionic lines are fitted with Gaussian profiles. 

We are mostly interested in the “residual” dust emission obtained by subtracting the 
PAH and ionic emission lines from the observed Spitzer/ IRS spectra. For illustration, we 
show in Figure [H the “residual” dust emission of IRAS F10398+1455, IRAS F21013-0739 
and SDSS J0808+3948 obtained by subtracting the PAH and ionic emission lines determined 
from the PAHFIT approach. We also show in Figured] the “residual” dust emission of three 
galaxies obtained with the spline approach. For F10398+1455, the overall “residual” dust 
emission profiles resulting from these two approaches are in close agreement with each other, 
except that the 9.7 /im silicate emission feature determined from the PAHFIT method is 
somewhat narrower than that from the spline method. In specific, the fluxes in the blue 
and red wings of the 9.7 pm. silicate emission feature resulting from the PAHFIT method 
are lower by <5% than that from the spline method. The flux differences between that de¬ 
rived from the spline method and that from PAHFIT are more pronounced in F21013-0739 


Table 1. Basic Parameters for the Three Spectroscopically Anomalous Galaxies IRAS F10398+1455, 

IRAS F21013-0739 and SDSS J0808+3948 


Sources 

R.A. 

Deck 

Redshift 

Classification 3, 

Program ID b 

M,° 

(Ms) 

Lm d 

(Lo) 

IRAS F10398 +1455 

10h42m33.32s 

+14d39m54.1s 

0.099 

AGN+SB 

40991 

io 10 - 7 

o 

p 

bi 

o 

IRAS F21013-0739 

21h03m58.75s 

-07d28m02.5s 

0.136 

AGN+SB 

40444 

io 10 - 9 

O 

p 

Cn 

CO 

SDSS J0808 + 3948 

08h08m44.27s 

39d48m52.36s 

0.091 

AGN+SB 

40444 

io 9 - 8 

IqIO.84 


a The optical classification was based on the Baldwin, Phillips & Terlevich (1981; BPT) diagram which plots 
the emission line ratio of [O III]/H/3 against [N Ilj/Ha (Y. Xie et ah 2015, in preparation). These three galaxies 
are classified as a composite of type 2 AGN and starburst. 

b The program ID from which the Spitzer/IRS spectra were taken. 

c Stellar mass (taken from the MPA-JHU SDSS DR7 catalog). 

d Infrared luminosity. 


Ox 
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Fig. 1.— The “residual” dust emission of IRAS F10398+1455 (a), IRAS F21013-0739 (b) and 
SDSS J0808+3948(c) obtained by subtracting from the Spitzer /IRS spectrum (black dotted line) the PAH 
and ionic emission lines determined from the spline approach (red solid line) or from the PAHFIT approach 
(blue dashed line). In the spline approach, the continuum underneath the PAH emission features is fitted 
with a spline function. In the PAHFIT approach, the continuum underneath the PAH emission features 
is approximated as a sum of starlight and blackbodies of different temperatures derived from the PAHFIT 
software of Smith et al. (2007). 
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and SDSS J0808+3948 than F10398+1455, with the spline-based 9.7//rri emission feature 
being broader and less prominent @ With dliiF u /dlnX & 2.7 and 2.6 for F21013-0739 and 
SDSS J0808+3948 respectively, the PAHFIT method results in a less steep ~ 5-8 //rri con¬ 
tinuum, but still much steeper than that of AGN (dlnF y /(ilnA~0.8). The relatively high 
flux density drawn from the spline approach could be due to an incomplete subtraction of 
the PAH wings. We will discuss in §5. II the differences in our modeling results derived from 
the two residual spectra. 


3. Dust Model 

To model the dust emission, we need to specify the dust chemical composition, size, 
and temperature. For dust composition, we will consider amorphous silicate and graphite or 
amorphous carbon, although other dust species (e.g., SiC, oxides) may be present in AGN 
tori (e.g., see Laor & Draine 1993, Markwick-Kemper et al. 2007, Koher & Li 2010). 

In the Milky Way diffuse ISM the dust grain sizes range from a few Angstroms to a few 
tenth micrometers (see Li 2004a). In AGN the dust is expected to be skewed toward a larger 
average grain size due to the preferential destruction of small dust in the harsh environments 
of AGN (see Li 2007), and probably also because of grain growth in the dense circumnuclear 
regions (see Maiolino et al. 2001). For simplicity, we will not consider size distributions, 
instead, we will only consider seven discrete sizes: a = 0.1, 0.5, 1.0, 1.5, 2.0, 5.0 and 10.0 /nn 
where a is the radius of the dust grains (assumed spherical). The dust sizes considered here 
are broad enough to cover both sub-/un-sized dust typical of the Galactic diffuse ISM (see Li 
2004a) and large, /im-sized dust typical of AGN (e.g., see Li et al. 2008, Smith et al. 2010, 
Lyu et al. 2014, Z. Shao et al. 2015, in preparation). 

The dust in AGN tori is expected to have a range of temperatures. In the inner “wall” 
region of the torus, the dust temperature reaches (>1500 K, the sublimation temperature of 
silicate and graphite materials. The dust temperature decreases when moving away from the 
inner “wall.” To fully account for the temperature distribution, we need to perform radiative 
transfer calculations for the dust in AGN tori heated by the central engine. For simplicity, 
we will only consider two temperatures — a warm component and a cold component — to 
represent the temperature distribution. 


2 We note that, although the spline-based 9.7 /irn silicate emission feature of SDSS J0808+3948 is almost 
too weak to be noticeable, the 18 ym O-Si O silicate bending feature is clearly seen (see Figure 2 of Xie et 
al. 2014). 
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We assume that the observed IR emission results from silicate dust and carbon dust 
(i.e., graphite or amorphous carbon). For each dust species, we consider a warm component 
with temperature T w and a cold component with temperature T c . Assuming the torus is 
optically thin in the IR, we model the dust IR emission as 

Fu x ^ ^ \^B p {T w i) x K a bs,?.(A) x Af W) j T B u (T C i ) x ACabs,*^) x M c , (1) 


where the sum is over the two dust species (silicate and graphite or amorphous carbon), 
d is the luminosity distance of the object, Kabs,*^) is the mass absorption coefficient or 
opacity (in unit of cm 2 g -1 ) of dust of type i, B U {T ) is the Planck function of temperature 
T at frequency u, T Wji and T C) , are respectively the temperatures of the warm and cold 
components of dust of type i, and and M c ^ are respectively the masses of the warm and 
cold components of dust of type i. For a given composition and size, we calculate the mass 
absorption coefficient K a b s (^) from Mie theory (Bohren & Huffman 1983) using the refractive 
index of the corresponding dust material. 

To summarize, there are in total eight parameters in fitting the observed IR emission: 
the temperature (T®) and mass (M®) for the warm silicate component, the temperature 
(I?) and mass (M®) for the cold silicate component, the temperature (T^) and mass (M°) 
for the warm carbon dust component, and the temperature (T c c ) and mass (M c c ) for the 
cold carbon dust component. We obtain the best fit for each galaxy using the MPFIT code, 
an IDL ^-minimization routine based on the Levenberg-Marquardt algorithm (Markwardt 
2009). 

We require the dust temperatures not to exceed the sublimation temperature (T sub i ~ 
1500 K) of silicate and graphite materials. For the dust mass, we could apply the elemental 
abundances to constrain the mass ratio of the silicate component to the carbon dust com¬ 
ponent. Let [X/H] tot be the total abundance of element X relative to H, [X/H] gas be the 
amount of X in the gas phase, and [X/H] d t be the amount of X contained in dust. As an 
element will be in the gas-phase and/or locked up in the solid-phase, one naturally obtains 
[X/H] dust = [X/H] tot — [X/H] gas . Let fix be the atomic weight of X (/ix~ 12, 16, 56, 24, and 
28 for carbon, oxygen, iron, magnesium, and silicon). If we assume a stoichiometric com¬ 
position of Mg 2 X Fe 2 (i_ x )Si 0 4 for silicate (i.e., each silicon atom corresponds to four oxygen 
atoms), the mass ratio of the carbon dust component to the silicate component is 


Afcarb/Af s il — 


_ ^c{[C/H] to ,-lC/H] SM } _ 

hFe [Fe/H] dust + /i Mg [Mg/H] dust + /igi [Si/H] dust + 4x/i 0 [Si/H] dust 


( 2 ) 
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If we assume that (i) these galaxies have the solar C, Fe, Mg and Si abundancesjj (ii) all Fe, 
Mg and Si elements are depleted in silicate dust, and (iii) ~ 37% of the C is in the gas-phase 
(i.e., [C/H] gas ~ 100ppm, Sofia et al. 2011), we derive M carb /M si i ~ 0.36. With C atoms all 
depleted in dust (i.e., [C/H]„. ag = Oppm), one obtains an upper limit of M carb /M si i ~ 0.57. A 
higher M carb /M si i ratio could be achieved if one assumes an enstatite (MgSiOa) composition 
for the silicate component: M carb /M si i « 0.59 for [C/H] gag = 100 ppm, or M carb /M si i ~ 0.95 
for [C/H] gas = Oppm. As we have no information about the elemental abundances and 
depletion for these three galaxies, in the modeling we will allow the mass ratio to vary 
between 0.2 < M carb /M si i < 2. 

We estimate the uncertainties for the eight model parameters by performing Monte- 
Carlo simulations. For the Spitzer/IRS spectrum of each source, we assume that the flux 
density error statistically follows a normal distribution. The dispersion is characterized by 
the observed la error, composed of the statistical and systematic errors. The latter arises 
from the flux differences between the two nods of the Spitzer /IRS spectra, the sky background 
contamination, and the Spitzer/IRS pointing and flux calibration errors (Leboutciller et al. 
2011)0 We generate a new “observational” spectrum through randomly sampling a point at 
each wavelength from the normal distribution. We then model the new spectrum and derive 
a set of model parameters. We conduct 100 simulations for each source as the parameters 
derived from 10,000 simulations only slightly differ from that derived from 100 simulations. 
The final model spectrum is calculated from the median values of the model parameters. 
The error of each parameter is derived from the standard deviation of 100 simulations. 

The mass absorption coefficient or opacity /c abs (zi) depends on the dust size, shape, and 
composition through the dielectric function or index of refraction of the dust material. For 
graphite, we take the index of refraction of Draine & Lee (1984). For amorphous carbon, 
we take the index of refraction of Rouleau & Martin (1991). For silicate, we consider a 
range of compositions: (i) “astronomical silicates” of Draine & Lee (1984), (ii) pyroxene 
Mg x Fei_ x Si 03 with x = 0.4, 0.7, 1.0 (Dorschner et al. 1995), and (iii) olivine Mg 2 x Fe 2 (i- x )Si 04 
with x = 0.4, 0.5 (Dorschner et al. 1995). In Figured we show the mass absorption coefficient 


3 [C/H] 0 « 269 ±31 ppm, [Fe/H] 0 « 31.6 ±2.9ppm, [Mg/H] 0 « 39.8±3.7ppm, and [Si/H] 0 « 32.4 ± 
2.2ppm (Asplund et al. 2009). 

4 We note that the signal-to-noise ratio (SNR) of the Spitzer /IRS spectrum in the ~ 5-14.5 //in wavelength 
interval is lower than that in the interval of ~ 14.5-38 /im. This was caused by the different observational 
modules (i.e., the 5-14.5 /im Short Low IRS module with a slit width of 376, and the 14.5-38 /an Long Low 
IRS module with a slit width of 1172). To fit the 5-8 //m continuum emission and the 9.7/im silicate emission 
feature, we arbitrarily increase the weights for the data points at 5-8 and 8-14.5 /im by ten and two times, 
respectively, compared to that at 14.5-38 /im. 
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ftabs(A) of each dust species of different sizes spanning the wavelength range of 2-40 pm. 

For all silicate species, the mass absorption coefficients K a bs(^) hi the 2-40 pm wave¬ 
length range are featured by a relatively smooth continuum at ~ 2-8 pm and two prominent 
bumps peaking around 9.7 and 18 pm. The ~ 2-8 pm continuum increases with the dust 
size. The 9.7 and 18 pm features also vary with the dust size: they substantially broaden 
and shift to longer wavelengths for a > 1.5 pm§ For amorphous carbon and graphite, the 
~ 2-8 pm continuum mass absorption coefficients flatten off with increasing dust size. 

The ~ 2-8 pm continuum mass absorption coefficients are sensitive to the silicate com¬ 
position. As illustrated in Figure [3}r, for a given dust size, the ~2-8 pm continuum opacity 
increases significantly with the iron fraction of the dust. In general, amorphous olivine has 
a higher ~ 2-8 pm continuum opacity than pyroxene. We note that the dielectric functions 
of the Draine & Lee (1984) “astronomical” silicate were synthesized to be “dirty” so that it 
is absorptive in the near-UV and optical wavelength ranges and also has a high ~ 2-8 pm 
continuum opacity. In Figure [3}o, we compare the 2-40 pm opacities of graphite and amor¬ 
phous carbon. In the ~ 5-8 pm wavelength range of interest here, amorphous carbon has a 
higher opacity than graphite. 


4. Results 

A first glance of the mass absorption coefficient /t a b s (A) profiles shown in Figures [21 [3] 
would lead one to speculate that the steep ~ 5-8 pm emission continuum observed in the 
three spectroscopically anomalous galaxies could be caused by silicate dust with an iron- 
poor pyroxene composition; and for a given composition, small, sub-pm-sized dust may be 
preferred over pm-sized dust (see Xie et al. 2014). In this section we apply the dust model 
elaborated in 1J3] to fit the observed Spitzer /IRS spectra of these galaxies by extensively 
exploring the parameter space of eight model parameters (T®, T c s , T^, T c c , M®, M®, M^, and 
M c c ) and different compositions (“astronomical silicate” [DL84], pyroxene MgSi0 3 [pyrO], 
pyroxene Mg 0 7 Feo. 3 Si0 3 [pyr3], pyroxene Mg 0 . 4 Fe 0 . 6 SiO 3 [pyr6], olivine MgFeSi0 4 [olv5], 
olivine Mgo.gFei^SiCL [olv6], amorphous carbon [ac], and graphite [gra]). We will focus on 
fitting the “observed” spectra obtained by subtracting from the Spitzer /IRS spectra the 
PAH and ionic emission lines determined with the PAHFIT method. The model fits to the 
“observed” spectra obtained with the spline method will be discussed in H5.ll 


5 For a < 0.5 (im (and even up to a < 1 /im), the 9.7 and 18 /im features are insensitive to the dust size 
since grains with a < 0.5 pm satisfy the Rayleigh scattering condition of 27ra/A <C 1 at these spectral bands 
(Bohren & Huffman 1983). 





2015.7.14.28: DRAFT 


11 






Fig. 2.— Mass absorption coefficients or opacities K a bs(A) for (a) “astronomical silicates” (Draine & Lee 
1984), (b) amorphous olivine Mgo. 8 Fei. 2 Si 04 (Dorschner et al. 1995), (c) amorphous olivine MgFeSiC >4 
(Dorschner et al. 1995), (d) glassy pyroxene MgSiC >3 (Dorschner et al. 1995), (e) glassy pyroxene 
Mgo. 7 Feo. 3 Si 03 (Dorschner et al. 1995), (f) glassy pyroxene Mgo. 4 Feo. 6 Si 03 (Dorschner et al. 1995), (g) 
amorphous carbon (Rouleau & Martin 1991), and (h) graphite (Draine & Lee 1984). All dust species are 
assumed spherical with a range of radii: a = 0.1/rm (red solid), a = 0.5^m (blue dashed), a=l/im (green 
dotted-dashed), and a = 2fim (purple dotted). 
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Fig. 3. — Top panel (a): comparison of the 2-40/mi opacity of the Draine & Lee (1984) “astronomical 
silicate” (“DL84,” solid red) with that of pyroxene MgSiC >3 (“pyrO,” dotted orange, Dorschner et al. 1995), 
pyroxene Mg 0 . 7 Feo. 3 Si 03 (“pyr3,” dashed purple, Dorschner et al. 1995), pyroxene Mgo. 4 Feo. 6 SiC >3 (“pyr6,” 
dotted-dashed green, Dorschner et al. 1995), olivine MgFeSiC >4 (“olv5,” dashed blue, Dorschner et al. 1995), 
and olivine Mgo.sFei. 2 Si 04 (“olv6,” dotted magenta, Dorschner et al. 1995). All silicate dust species are 
taken to be spherical. All dust species are assumed spherical with a radius of a = 0.1 /im. Bottom panel 
(b): comparison of the 2-40/mr opacities of graphite (Draine & Lee 1984) of a = 0.1 /jm (solid blue) and 
a = 1 /im (dashed blue) with that of amorphous carbon (Rouleau & Martin 1991) of a = 0.1 /jm (solid red) 
and o = l /mi (dotted-dashed red). 










Table 2. Model Parameters for Three Spectroscopically Anomalous Galaxies (IRAS F10398+1455, IRAS F21013-0739, and 
SDSS J0808+3948) and Two Comparison Sources (Quasar Average Spectrum, and ULIRG IRAS FSC 10214+4724) 


Object 

Composition 

Dust Size 

rjiS 

rjp C 

rpS 

rp C 

(10 3f g) 

M S 

(nd^g) 

M s 

(10 3f g) 

M c 

(10 3? g) 

X 2 /dof 



(H 

(K) 

(K) 

(K) 

(K) 


F10398+1455“ 

DL84 + ac 

0.1 

271+1.6 

389+2.8 

92+0.2 

117+0.1 

14.0+0.5 

3.1+0.2 

5.9+0.03 

11.9+0.07 

38 

F10398+1455 6 

DL84 + ac 

0.1 

248+4.9 

295+1.3 

84+0.5 

115+0.3 

15.4+2.1 

27.0+4.3 

7.4+0.08 

14.8+0.2 

42 

F10398+1455“ 

olv6 + ac 

0.1 

303+2.0 

323+2.2 

76+0.3 

113+0.1 

8.8+0.3 

10.9+0.5 

8.7+0.04 

17.4+0.09 

25 

F10398+1455“ 

pyr3 + ac 

1.0 

248+3.0 

259+1.5 

74+0.8 

109+0.3 

16.0+1.1 

20.6+1.9 

9.1+0.2 

18.3+0.4 

24 

F10398+1455“ 

DL84 + gra 

0.1 

291+4.2 

319+2.6 

120+0.6 

141+0.7 

7.6+0.5 

15.2+1.1 

1.7+0.05 

2.7+0.1 

36 

F10398+1455“ 

DL84 + ac 

1.0 

350+30.0 

284+3.9 

93+0.4 

113+1.1 

3.7+2.9 

7.4+6.0 

6.1+0.1 

12.2+0.2 

59 

F10398+1455“ 

DL84 + gra 

1.0 

395+13.8 

385+5.7 

112+0.1 

121+0.4 

2.4+0.3 

1.1+0.3 

1.8+0.02 

3.1+0.08 

30 

F10398+1455“ 

olv6 + ac 

1.0 

346+4.0 

285+1.6 

81+0.4 

111+0.3 

4.5+0.2 

9.0+0.4 

8.0+0.1 

16.0+0.2 

30 

F10398+1455“ 

DL84 + ac 

0.1 

244+3.3 

357+3.2 

92+1.0 

115+0 

23.7+2.1 

5.4+0.6 

6.2+0.2 

12.5+0.3 

42 

F10398+1455“ 

DL84 + ac 

1.5 

410+4.5 

374+4.1 

93+0.2 

109+0.1 

2.4+0.1 

0.6+0.05 

6.2+0.03 

12.4+0.05 

99 

F10398+1455“ 

olv6 + ac 

1.5 

358+13.4 

354+4.5 

84+1.1 

108+0.2 

4.0+0.5 

1.3+0.3 

7.6+0.3 

15.0+0.5 

51 

F10398+1455“ 

pyr3 + ac 

1.5 

331+3.5 

385+3.0 

78+1.6 

109+0.4 

4.5+0.2 

0.9+0.04 

8.3+0.3 

16.6+0.7 

27 

F10398+1455“ 

DL84 + ac 

2.0 

419+3.1 

334+3.0 

99+0.1 

99+0.1 

2.8+0.06 

1.2+0.1 

6.0+0.02 

12.0+0.04 

201 

F10398+1455“ 

olv6 + ac 

2.0 

464+50.0 

324+14.0 

89+0.2 

99+0.5 

1.5+1.3 

2.6+2.4 

7.9+0.1 

15.8+0.2 

168 

F10398+1455“ 

pyr3 + ac 

2.0 

641+8.4 

382+2.4 

89+0.1 

102+0.1 

0.6+0.02 

1.1+0.04 

7.5+0.02 

15.0+0.05 

101 

F21013-0731“ 

DL84 + ac 

0.1 

249+4.0 

361+5.0 

61+0.4 

110+0.5 

24.6+2.1 

4.9+0.4 

15.0+0.8 

19.8+1.8 

18 

F21013-0731 6 

DL84 + ac 

0.1 

240+0.0 

301+1.3 

74+0.8 

120+3.6 

16.2+0.5 

32.5+1.0 

20.7+4.4 

10.3+4.0 

21 

F21013-073l“ 

olv6 + ac 

0.1 

334+5.0 

318+4.6 

58+0.4 

113+0.6 

6.3+0.4 

12.0+0.1 

49.0+3.9 

16.7+2.5 

17 

F21013-0731“ 

pyr3 + ac 

1.0 

251+9.2 

255+3.2 

66+3.7 

110+1.1 

16.4+4.5 

21.7+6.9 

26.0+9.4 

16.3+19.3 

17 

F21013-0731“ 

DL84 + gra 

0.1 

267+0.5 

327+5.5 

111+0.4 

144+0.6 

14.4+0.02 

12.2+1.4 

1.7+0.003 

2.6+0.03 

19 

F21013-073l“ 

DL84 + ac 

0.1 

270+4.3 

374+5.3 

74+2.0 

115+0.0 

15.0+1.4 

3.7+0.6 

9.7+1.1 

15.0+2.6 

20 

J0808+3948“ 

DL84 + ac 

0.1 

265+55.0 

363+50.0 

84+11.7 

110+2.0 

1.7+0.9 

0.3+0.2 

1.7+0.9 

2.4+1.4 

0.3 

J0808+3948 6 

DL84 + ac 

0.1 

200+0.0 

285+3.5 

78+11.6 

106+2.0 

3.2+0.2 

6.3+0.5 

2.1+0.4 

3.0+0.3 

0.7 

J0808+3948“ 

olv6 + ac 

0.1 

297+58 

383+53.0 

65+5.9 

108+3.1 

1.3+1.0 

0.3+0.3 

2.3+1.4 

3.7+2.6 

0.2 

J0808+3948“ 

pyr3 + ac 

1.0 

239+24.5 

294+26.0 

67+5.9 

105+2.4 

2.5+1.5 

0.5+0.4 

2.1+0.7 

4.1 + 1.5 

0.8 

J0808+3948“ 

DL84 + gra 

0.1 

265+56.0 

379+66.0 

106+6.5 

129+7.4 

1.7+1.3 

0.3+0.3 

0.5+0.2 

0.5+0.3 

0.5 

J0808+3948“ 

DL84 + ac 

0.1 

285+34.0 

384+35.0 

87+3.6 

115+0.0 

1.1+0.5 

0.2+0.1 

1.7+0.3 

1.8+0.4 

0.3 

Quasar 

DL84 + gra 

1.5 

350+32.0 

644+5.8 

40+6.0 

149+1.3 

- 

- 

- 

- 

0.42 

FSC 10214+4724 

DL84 + gra 

0.5 

394+78 

377+132 

72+35 

138+11 

5.2 x 10 a 

2.4 x 10** 

1.9 X 10* 

5.6 X IQ- 1 

5.8 


a Model parameters derived from the PAHFIT-based IR emission, 
kModel parameters derived from the spline-based IR emission. 
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Fig. 4. — Comparison of the model emission (solid green line) of IRAS F10398+1455 (a), IRAS F21013-0739 
(b), and SDSS J0808+3948 (c) with the “observed” spectra (black open squares) obtained by subtracting 
from the Spitzer /IRS spectra the PAH and ionic emission lines determined with the PAHFIT method. The 
model consists of four components: warm silicate (red dashed line), cold silicate (magenta dot-long dashed 
line), warm amorphous carbon (blue dotted-dashed line), and cold amorphous carbon (orange long dashed 
line). All dust components are taken to be 0.1 fim in radii. The DL84 “astronomical” silicate is adopted for 
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4.1. Dust Composition 

We first consider the DL84 “astronomical” silicate and amorphous carbon, and take all 
dust components to have the same grain size of 0.1 /im in radii. The effects of dust size 
will be discussed in 1 14.21 In Figure 0] we show the model fits to the “observed” spectra 
of IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948. The overall fits are 
very good except for SDSS J0808+3948; the model 9.7 /irn silicate profile is a little bit too 
broad and its peak wavelength is a little bit too short. Also, for IRAS F10398+1455 and 
IRAS F21013-0739, the model slightly overpredicts the red wing of the 9.7 /irn feature. For all 
three sources, the ~ 5-8 fi m emission continuum is dominated by warm amorphous carbon, 
and the emission at A > 14 /irn is dominated by cold amorphous carbon. 

We then consider models consisting of amorphous olivine and amorphous carbon. As¬ 
suming a size of a = 0.1 fi m for all dust components, we find that Mgo.gFeoSiCR (together 
with amorphous carbon) provides an excellent fit to all three sources (see Figure [5]). Sim¬ 
ilarly, we have also considered models consisting of amorphous pyroxene and amorphous 
carbon. As shown in Figure [6], excellent fits are achieved with Mgo. 7 Feo. 3 Si 03 of a = 1.0 /mi. 
This is because the 9.7/irn feature of pyroxene dust is appreciably narrower than that of 
amorphous olivine and “astronomical” silicate (see Figure [3]). Therefore, models consisting 
of a = 0.1 /irn pyroxene and amorphous carbon produce too narrow a 9.7/irn feature. On 
the other hand, if the dust is too large (say, a = 2.0/mi), the model 9.7/irn feature will be 
too broad to be consistent with the observed spectra (see 1 14.21) . 

We have also considered graphite. As shown in Figure [71 models consisting of graphite 
and “astronomical silicate” of a = 0.1 /irn closely fit the observed IR emission of all three 
galaxies, with the fitting y 2 at a similar level as that of the “astronomical silicate” /amorphous 
carbon model (see Figure 0]). 


4.2. Dust Sizes 

Xie et al. (2014) speculated that the steep ~5-8 /im emission continuum seen in these 
spectroscopically anomalous galaxies could arise from sub-/nn-sized silicate dust instead of 
/im-sized silicate dust since the latter displays a larger opacity in the ~5-8 /im wavelength 
range (see Figure [2]). However, as demonstrated in 1 14.11 the ~5-8/mi emission continuum 
is dominated by the carbon dust component. Therefore, the ~5-8 /mi continuum is not a 
powerful constraint on the silicate dust grain size. Instead, the 9.7/irn feature allows us 
to place a stronger constraint. Indeed, as shown in Figures I8HTU1 for models consisting of 
different silicate and carbon dust species (“astronomical silicate,” olivine, pyroxene, graphite, 
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Fig. 5. — Same as Figure 0] but with the DL84 “astronomical silicate” replaced by amorphous olivine 
Mg 0 . 8 Fei. 2 SiO 4 (“olv6”). 
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Fig. 6.— Same as Figure [I] but with the DL84 “astronomical silicate” replaced by amorphous pyroxene 
Mgo. 7 Feo. 3 Si 03 (“pyr3”) and all dust components taken to have a size of a = 1.0/im. 
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Fig. 7. — Same as Figure@]but with the amorphous carbon component replaced by graphite. 
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Fig. 8. — Upper panel (a): Same as Figure [I| a) but with a = 1.0/rm for all dust components. Middle 
panel (b): Same as Figure EJa) but with a = 1.0/im for all dust components. Bottom panel (c): Same as 
Figure [5ja) but with a = 1.0 /jm for all dust components. 
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Fig. 9. — Upper panel (a): same as Figure @|a) but with a = 1.5 /im for all dust components. Middle 
panel (b): same as Figure [SJa) but with a = 1.5 /im for all dust components. Bottom panel (c): same as 
Figure [Gja) but with a = 1.5 /im for all dust components. 
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Fig. 10. — Same as Figure|H]but with a = 2.0 fim for all dust components. 
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and amorphous carbon), the observed IR emission could still be closely fitted with a = 
1.0 /mi, but not with a = 1.5 /mi or a = 2.0/mr. This is easy to understand: as illustrated 
in Figure [21 the 9.7 //m profile is insensitive to dust size for a < 1 /irn, but it considerably 
broadens for larger sizes. The best-fit models with a = 1.0 /irn lead to an appreciably higher 
temperature for the warm silicate component (T®) compared with that of the best-fit models 
with a = 0.1 jum. Consequently, the mass required by the warm silicate component (M^) 
drops by a factor of ~ 3-5. This can be understood in the sense that, compared with silicate 
dust of a = 0.1 /rm, the “9.7 /i m” silicate feature of the opacity K a bs( z/ ) profile of silicate 
dust of a = 1.0/mi slightly shifts to a longer wavelength. To compensate this “redshift,” 
one requires an increased temperature which blueshifts the peak wavelength of the Planck 
function B U (T), and therefore the resulting emission spectrum which is the product of B U (T ) 
and K a bs(r) remains at the observed peak wavelength. The effects of grain size on the derived 
dust temperature are manifested in Figures [91 and [TUI 


5. Discussion 

We have shown in 1J4] that a simple model consisting of four dust components (warm 
and cold silicate dust, and warm and cold carbon dust) closely reproduces the observed 
IR emission of all three spectroscopically anomalous galaxies, including the steeply rising 
~5-8 /mi emission continuum and the 9.7 /mi emission feature. The model-fitting results are 
insensitive to the exact silicate or carbonaceous dust composition. As shown in Figures [4HZ1 
models consisting of “astronomical” silicate, amorphous olivine, or pyroxene, combined with 
amorphous carbon or graphite, are all capable of successfully fitting the observed IR emission. 
The ~5-8 /mi emission continuum is dominated by the warm carbon dust component with 
T in the range of ~250-400K. The 9.7/irn silicate emission is predominantly due to silicate 
dust of T ~ 200-400 K. The emission at A > 14 /mi is dominated by the cold carbon dust 
component mostly with T ~ 110-120 K. The 9.7 /mi emission feature constrains the silicate 
dust size to not exceed ~1.0/mi (i.e., a < 1.0 /mi, see Figures I8HTUD. The mass ratio of the 
warm carbon dust to the warm silicate dust ranges from ~0.2 to ~ 2.0, with a 

mean ratio of ~0.98. Similarly, the mass ratio of the cold carbon dust to the cold silicate 
dust (M®/M c c ) ranges from ~0.34 to ~2.0, with a mean ratio of ~ 1.5. The total dust 
mass is dominated by the cold components, with the warm components only accounting for 
< 0.15% of the total dust mass. 

To investigate the significance of the dust temperature and mass yielded from our dust 
IR emission modeling, for all three sources we fix the temperature of the cold carbon dust 
component to be T c c = 115 K and then probe the dust temperature and mass variations for 
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the other components. Figure [TT] shows the model spectra for all three sources obtained from 
the mixture of “astronomical silicate” and amorphous carbon (“DL84 + ac”) with a = 0.1 /irn 
and T c c = 115 K. The model spectra are in close agreement with the observed spectra. 
Compared with the same model but with T c c treated as a free parameter (see Figure 2]), the 
temperatures derived for the other three components (T®, T^, and T c s ) vary within < 10%, 
while the derived masses (M®, M^, and M®) vary within < 75%. 

To test the reasonableness of the dust mass derived in §H] f° r each source we compare 
the dust mass with the stellar mass. On average, the dust-to-stellar mass ratio of these three 
galaxies is ~ 10~ 5 , much smaller than that of the Milky Way (~ 10 -3 ; see Li 2004b). This 
ratio appears reasonable since the mid-IR emission considered here only probes the dust of 
T > 100 K (see Table [2]), while the bulk dust mass is in dust with T < 100 K in the host 
galaxy which emits in the far-IR and escapes from detection by Spitzer/ IRS. 


5.1. PAHFIT Versus Spline 

The results presented in §|4] are derived by fitting the “observed” 1R emission obtained 
by subtracting from the Spitzer/IRS spectra the PAH and ionic emission lines determined 
with the PAHFIT method. As illustrated in Figured] the 1R emission spectra derived from 
the spline method show a relatively broad and weak 9.7 //m silicate emission feature in com¬ 
parison with that from PAHFIT. Furthermore, the spline method yields a steeper slope for 
the ~ 5-8/mi continuum than PAHFIT (see 1|2]). We have also modeled the spline-based 
“observed” 1R emission of all three spectroscopically anomalous galaxies. As illustrated 
in Figure d21 the simple four-component model consisting of warm/cold silicate dust and 
warm/cold carbon dust closely reproduces the “observed” 1R emission obtained from the 
spline method. The temperatures and masses of the cold components for all sources do 
not change much (<10% for T, and <30% for M; see Table [2]) since the PAHFIT- and 
spline-based IR emission spectra exhibit little difference beyond A ~ 14 //m where the cold 
components dominate the radiation. In contrast, the temperature of the warm silicate com¬ 
ponent derived from the spline method drops by ~ 10-25% and its mass changes by a factor 
of ~1.1-1.9. For the warm carbon component, on average, its temperature drops by ~20%^| 
and the corresponding mass increases by a factor of ~ 8.7, ~6.6 and ~21 for F10398+1455, 
F21013-0739 and SDSS J0808+3948, respectively. 


6 Compared with PAHFIT, the spline method yields a steeper and higher ~5-8/xm continuum (see Fig¬ 
ure |T| and therefore a lower temperature and larger mass for the warm carbon component which dominates 
the ~ 5-8 /im continuum. 
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Fig. 11. — Same as Figure @] but with the temperature of the cold carbon component (orange long dashed 
line) fixed at T c c = 115 K (see !j5]). 
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Fig. 12. — Same as Figure[]]but with the “observed” spectra (black open squares) obtained by subtracting 
from the Spitzer /IRS spectra the PAH and ionic emission lines determined with the spline method (see (JT|). 
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Fig. 13.— Comparison of the observed IR emission of IRAS FSC10214+4724 (black open squares; Teplitz 
et al. 2006) with the model emission (solid green line) which is a sum of four components: warm DL84 silicate 
(red dashed line), cold DL84 silicate (magenta dot-long dashed line), warm graphite (blue dotted-dashed 
line), and cold graphite (orange long dashed line). All dust components have a size of a = 0.5 fj, m. 


5.2. Comparison with ULIRG and Quasars 

The anomalous IR emission (i.e., a steep ~ 5-8 pm emission continuum combined with 
silicate emission) seen in IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948 
is rarely seen in other sources. To the best of our knowledge, the only similar source is 
IRAS FSC 10214+4724, a lensed, starburst-powered ULIRG at a redshift of z ~ 2.29 (Teplitz 
et al. 2006). The major difference between IRAS FSC 10214+4724 and our three galaxies 
is that IRAS FSC 10214+4724 lacks the PAH emission features which are prominent in our 
three galaxies. We have also modeled the IR emission of IRAS FSC 10214+4724. As shown 
in Figure fTHl a mixture of warm/cold “astronomical silicate” and warm/cold graphite of 
a ~ 0.5pm closely fits the observed IR emission. The derived parameters (see Tabled are 
similar to that of our three galaxies. 
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Fig. 14.— Same as Figure H3l but for the average spectrum of quasars (see Hao et al. 2007). 


We have also modeled the average spectrum of quasars (Hao et al. 2007)13 As shown in 
Figure [14l the average quasar spectrum exhibits a flat ~5-8 /irn emission continuum and a 
broad silicate emission feature at 9.7//m. The best-fit is achieved with dust of a ~ 1.5 /irn, 
with T® ~ 350 K for warm silicate, ss 644 K for warm graphite, T c s ^ 40 K for cold 
silicate, and T c c « 149 K for cold graphite. The flat ~ 5-8 /irn emission continuum seen in 
quasars is caused by warm graphite which is about twice as hot as that in our three galaxies. 

In principal, it would be desireable to model the IR emission spectra of starbursts to 
which the steep ~ 5-8 /im continuum emission and the PAH emission features of our three 
galaxies bear a close resemblance. However, to obtain any meaningful information about the 
properties of the dust (e.g., size, composition, temperature) from modeling the IR spectra 
of starbursts, one has to make assumptions concerning the distribution of stars and dust 
within the starburst. We therefore intend not to model the IR spectra of starbursts as we 
know little about how stars and dust are spatially distributed. 


7 The averaged spectrum of quasars is constructed with individual quasar spectra normalized at A = 
14.5/rm, in order to avoid the most luminous quasar dominating the resulting spectrum. 
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Y. Xie et al. (2015, in preparation) examine the multi-wavelength data of these three 
galaxies to investigate their powering sources. They find that F21013-0739 and SDSS J0808+3948 
have detections in hard X-rays (Jia et al. 2011) which indicates the presence of AGN. How¬ 
ever, the diagnostics from the UV, optical, mid-IR, and radio do not show any strong evi¬ 
dence of AGN. They conclude that AGN may be present in these galaxies, but they are too 
young to be detected in all wavelength bands. The lower dust temperature (T ~ 250-400 K) 
and smaller grain size (with an upper limit of a < 1.0 /im) inferred for our three galaxies 
compared with that of quasars could be due to the fact that our galaxies may harbor a 
young/weak AGN which is not maturely developed yet. Compared with quasars, such a 
young AGN emits much less and softer in the X-ray and UV. Therefore, sub-/an-sized grains 
could survive against sublimation and photo-sputtering in our galaxies while only /im-sized 
grains could survive in quasars. Exposed to a young AGN with a lower UV luminosity Luv; 
the dust temperature is expected to be lower since T oc L^y i+l3 \ where ft is the dust opacity 
power-law exponent in the far-IR (i.e., K a b s oc A -/3 ). Indeed, the bolometric luminosity of 
our galaxies is in the order of ~ 10 44 ergs~ 1 (Y. Xie et al. 2015, in preparation), about two 
orders of magnitude lower than that of typical quasars. 


6. Summary 

We have modeled the Spitzer /IRS spectra of three spectroscopically anomalous galaxies 
(IRAS F10398+1455, IRAS F21013-0739 and SDSS J0808+3948) in terms of a simple model 
consisting of a mixture of warm and cold silicate dust, and warm and cold carbon dust. The 
IR spectral characteristics of these galaxies are unique in the sense that they show silicate 
emission which is characteristic of AGN while they also show a steep ~ 5-8/im emission 
continuum and strong PAH emission features which are characteristic of starburst galaxies. 
In contrast, AGN exhibit a flat emission continuum at ~ 5-8 /an and lack the PAH emission 
features. The 9.7 and 18 /an silicate features seen in emission in these three galaxies are 
seen in absorption in starbursts. We find that our model with four components (cold and 
warm carbons, cold and warm silicates) closely reproduces the observed IR emission of all 
three galaxies, including the steep ~ 5-8/im emission continuum and the 9.7 /an silicate 
emission feature. The model-fitting results are insensitive to the exact silicate and carbon 
dust composition (e.g., “astronomical” silicate or laboratory analogs, olivine or pyroxene, 
iron-poor silicate or iron-rich silicate, amorphous carbon or graphite) and dust size, provided 
the dust size does not exceed one micrometer (i.e., a < 1.0/an). It is also found that the 
dust temperature is the primary cause in regulating the steep ~ 5-8 /an emission continuum 
and the prominent silicate emission feature. Compared with the average spectrum of quasars 
which exhibits a flat ~ 5-8 /an continuum, the steep ~ 5-8 /an continuum seen in our three 
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galaxies predominantly arises from warm carbon dust of temperatures of T ~ 250-400 K, 
much lower than that of quasars (T ~ 640 K). The 9.7 /mi silicate emission of the quasar 
average spectrum is best fitted with silicate dust of size a ~ 1.5 /mi, larger than the upper 
size limit of a < 1.0 /mi inferred for our three galaxies. The lower dust temperature and 
smaller dust size inferred for our three galaxies compared with that of quasars could be 
due to the fact that our galaxies may harbor a young/weak AGN which is not maturely 
developed yet. To test the simple dust model used in current work, we would like to present 
a series of physical AGN models (e.g Clumpy models) for these three galaxies as well as the 
quasar average spectrum and ULIRG FSC 10214+4724 using the model parameters derived 
in this paper (R. Nikutta et ah 2015, in preparation). 
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